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PART (1) 


BOUNDARY LAYER IN COMPRESSIBLE FLUIDS 


The solution of flow problems in which the density 
is variable is in general very difficult; hence, every case 
in which an exact or even an approximate solution of the 
equations of the motion of compressible fluids can be obtained 
has considerable theoretical interest, Several authors noticed 
that the theory of the laminar boundary layer can be extended: 
to the case of compressible fluids moving with arbitrarily high 
velocities without encountering insurmountable mathematical 
difficulties. Busemam (Ref. 1) established the equations and 
calculated the velocity profile for one speed ratio. (By speed 
ratio is understood the ratio of the eirspeed to the velocity 
of sound.) Frankl (Ref. 2) also made an analysis of the same 
problem, however, it is complicated and depends on several 
arbitrary approximations. Von Karman (Ref. 3) obtained a 
first epproximetion by a simple but apparently not sufficiently 
exact calculation. Hence, in Section (I), e better method for 
the solution of the problem is developed. 

The houndary layar theory for very high velocities 
is not without practical interest. First, the statement can 
be found often in technical and semi-technical literature on 
rockets and similar high-speed devices that the skin friction 


becomes more and more significant at high speeds. Of course, 


“ls 


it is krown that with increasing Reynolds Number, the skin 
friction coefficient is decreasing, i.e., the skin friction 
becomes relatively small in denver iden with the drag produced 
by wave formation or direct shock. Since high-speed flight 
will be performed mostly at high eltitude where the air is of 
wees low density, so that the kinematic viscosity is large, 
the resulting Reynolds Number will be relatively smell in 
spite of the high speed. 

Another interestins point in the theory of the 
boundary layer in compressible fluida is the thernodynamic 
aspect of the problem. In the case of low speeds the in- 
fluence of the heat produced in the boundary layer can be 
neglected both in the calculation of tho drag and of the heat 
transferred to the wall. In the case of high speeds, however, 
the heat produced in the boundary layer is not negligible, 
but determines the direclion of heal Pluw. In Section (II) 

a few simple examples of heat flow through the boundary 
layer are discussed, 

it nes been found necessary in most parts of this 
analysis to make the assumption of laminar flow. “his assump- 
tion was found necessary because of the lamentable state of 
knowledge concerning the laws of turbulent flow of compressible 
fluids at high speeds. This assumption is somewhat justified 


by the fact that - as mentioned above - in many problems where 


=D a 


the ts of this paper can be applied, the Reynolds Number 
is relatively small, so that a cómmidonable portion of the 
boundary layer is probably, de facto, laminar. Ackeret (Ref .l) 
called ettention to the possibility that the stability 
conditions in supersonic flow might be quite different from 
those occuring in flow with low velocities. 

Recently Kichemann. (Ref.5) studied the stability 
of a linear profile near a wall under small sinuous disturbance, 
and shown that as the velocity ratio increases, the flow is 
unstable at increasing wave length of disturbance. However, he 
assumed that the gas is non=viscous and the er velocity is 
constant., Both assumptions tend to limit the usefulness of the 
theory, especially the later ore, Because the Constancy of 
velocity of gas implies the constancy of gas temperature, which 
anfortunately is far from the truth for a perfect gas as will 
be shown in later calculations of this section. In spite of 
these enoescninties: some calculations of this section, as 
will be pointed out, are also applicable to turbulent flow. 
In other cases, as in the calculation of drag, the assumption 
of laminar flow surely gives at least the lower limit of its 
value. 

Section (1) 
If the X -axis is taken along the plate in the 


directiou of Lhe free slream, Lhe Y eyxis perpendicular Lo the 


-j= 


plate, and & and Y indicate the MX and Y components 
of the velocity at any point, then the simplified equation of 


P š 7 i a F ? ee 
motion in the boundary layer is , te. ios oy” 


su J Pio = 4 (up) (1.1) 


where both the density $ and the viscosity /4 are variables. 


The equation of continuity in this case is 
2 2 = 
A (Pu) + y (Pv) =0. (1,2) 


A third equation determines the energy balance between 
the heat produced by viscous dissipation and the heat transferred 
by conduction and convection. With the same simplification as 


used in Eyse(1.1 and (1,2), one can write 


2 2 E 
ECT) ter Gl OS) E) 
(1.3) 
where Cp is the specific heet at constant pressure, and À 
is the coefficient of heat conduction. If Prendtl's number, 
bm JA is assumed to be equal to 1, then it can be 
easily shown that both Eqs. (1.1) and (1.3) can be satisfied by 
equating the temperature T to a certain parabolic function of 


the velocity 4 only. Indeed, introducing G7 = Hu) 


ale 


into equation (1,3) and replecing A by E Me » one 


obtains 


(se do AUGE) = Sql pe H tia pf P rare 


$ 
Hence Eq. (1.3) is reduced to Eq. (1.1) if A ak -7 


or 

where É, E G are constants, Denoting the wall 
temperature [u = 0) by Ly and remembering that 
T= ya for 4=(/ where (f = free 
stream velocity, C, C; can be expressed in terms 


+ 


of tg: Ta and 


Lifferentiating Eo. (1,4) om obtains 
FW) = Ulm SE) (1.5) 
ha I ULZ FA LAG Lor 


where the subscript wr refers to conditions existing at the 
Le 
surface of the plate. Now (34/24), is elways 


positive; therefore, if MDA > T -J 


heat is trensferred from the fluid to the wall; if 


[4-02] = u/n -1 
there is no heat transfer between the fluid and the well; and 


if 


/ Ad .. — f E 
[RYL IT < Turf Ta > 1 
heat is transferred from the wall to the fluid. If there is 


ro heat transfer, the energy content per unit mass 


(url ) + G Fa is constant in the whole 


rerion of the boundary layer. (Ref.6) 


The pressure being constant the relation between 


Pf and EA is, 


f = Soa (1.6) 


bm 


The expression for the viscosity based on the kinetic theory 


of gases is 


pe 
PERRA (1.7) 


However, the following formula is in closer agreement with 


experimental data 
7 07 
Ms Ma ( z/ (1.78) 


Busemann (Ref. 1) calculated the limiting case for 
x 2 
which [py JM = ( BE e 5 using Eq. 
(1.7) and found that for a high Mach's number, the velocity 
profile is approximately linear. Von Karman (Ref. 3), using 
the linser velocity profile, the integral relation between 


the friction end the momentum, and Eq. (1.7) found that 


A Frictional drag per unit width of plate 
fi (£U ys ) x length of plate 


z po Ll pq? 
0 LUX fr A J 


The dimensions quantity (O) shown in Table (1.1) is a 


+ 
4 (1.8) 


function of Mach's mmber only. 


However, if Eq. (1.78) is used, then 


Be fo 


-4.12 


«ob tim 
, O fe fst EM 


(1.88) 


Table 1.1 
M a 1 2 5 10 eo 
© 1.16 1.20 1.25 1.39 1.50 1.57 


It is evident that this linear approxmiation is not 
satisfactory for smell values of Mach's number. For M =0 , 
the case is the same as the Blasius solution (Ref. 7) for in- 
compressible fluids for which © is 1.328. 

To solve the problem more rigorously, one has to 
resort to Eqs. (1.1) and (1.2). By introducing the streem 


function y which is defined by 


E oyak 
ae 7 Y OX 


the equation of contimity, Eq. (1.2), is satisfied automatically. 
Now if in Eq. (1.1) ye is introduned as the independent 
variable as was done by von Mises (Ref. 8) in his simplification 


of the boundary leyer equation for incompressible fluids then 
há 2 
y LED it in UP 
Y g W, 
where »L is a coordinate measured in normal direction to 
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Y coordinate. Using these new coordinates, the following 
relations exist 


su Ta UG Wf 


ag iy “MELA, wp 


therefor Eq. (1.1) oan be written as 


es a LM A 


(1.90) 


Eq. (1.9a) can be put into non-dimensional form by 
introducing the following set of new quantities: 


uy 
= nfl, 


(1.9b). 
e (PIILU = lM 
e* = sg 


po = pf, 


where L is a convenient length, say length of the plate, 


is the corresponding Reynold's Number, then Eq. (1.9a) 
becomes 


f * out 
Tit = ya (Use a 


(1.9) 
gde: 


Eq. (1,9) can be further simplified by introducing 


a new dependent variable e = y* y nt » then 


“Ísis A (ut stut GE) 


This ean be a by the method of successive 


(1.10) 


approximations. As st and pr are functions of temperature 
only as shown in Egs. (1.6 and (1,7) or (1.78) and the temper- 
ners is a function of ut then by starting with the known 
Blasius' solution (Ref. 7) the right-hand side of Eq. (1.10) 


can be expressed in terms of E « Therefore, one can write 


ub tut HG) and Ege (1.10) bocomes 5 ge LHe] 


Cuusequently, the solution of Bge (1.10) is 


g 
F _f SH 
CIF o A 


where = e 


(1.11) 
and C is determined by the boundary condition that at Ç = 00 


uol + = Cs AC (lella) 


In the actual computation, twa methods of evaluating 
the integrals in Eqs. (1.11) and (l.lla) are used. For small 
A+ 
values of Ç ; €<02 , the function 4 and 
118) ere expended in a power series of S - Due to 


the uniform convergence of the power series for sufficiently 


= 10 


small values or € , the integration is carried out term by 
term. For values of Ç >02  , numerical integration is 
used. | 

À second approximation can be made based upon the 
value of u* obtained from Eg. (1.11). It has been found 
in the cases investigated that the third end fourth approxima- 
tion gives sufficient accuracy, if the velocity profile of next 
smaller Mach's mber is used as the starting point of calculat- 
ing the velocity profile of next larger Mach's number, 

Having computed the final 4 Pe Y 
corresponding to y* can be calculated in the following 
ways 


It is known from Eq. (1.9b) that 


oe” E 


Then remembering the definition of Wp » one has 


ESA 
ay? Anat 


AA 


However due to the smell slope of stream lines, 


Es 
y* dp 


=il< 


E | 
dy? = ss HE (1.128) 


NR y — / É Ln 
nt te X feo [ae (1.12) 
US, 


Here the expansion of et y* in a power 
series for small velues of a is especially useful due to the 
singularity of integrand at e =) . 


The skin friction oan be computed by using momentum 


Drag = /)= / E 4 ul U- u) df, 


Using Eqs (1.120), on hao 


low, 1.004 


AF La 
a di Eee ue 


thus 


De fis Ef fra 


Atk 
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But y A = EA at therefore 
i 6 
Hence 


Thus the 


skin friction coefficient can be computed as 
/ D 2/ (1- Wd 
= OOO ma “O 
= 7 = 
f FYL JR (1.13) 


The velocity profile, the temperature distribution, 
and the frietional drag coefficient are ealenleted for different 


values of the Mach's number of the free stream, for the case 
2 
[ve J11?=(Ta/T)- 1 


using the approximate viscosity relation of Eq. (1.78). The 
results are shown in Figs. 1.2 and 1.3. The velocity profiles 
for high speeds are very nearly linear, but it can be seen that 
the wall temperature for greater Mach's numbers is very high. 
If the free stream temperature is LOOF., then the wall tempera- 
ture will be 1600ºF., 3620F., 6540°F., and 10,170ºF. for 


-13- 


Mach's numboro of l, 6, B, end 10, REP Poe Therofore, 
there is no doubt that the law of viscosity as expressed by 
Eq. (1.78) will not hold. Algo. et such high temperatures, the 
heat transfer due to radiation cannot be neglected, The effect 
of radiation will be the equalization of gas temperature. In 
the extreme case of complete equalization, the temperature 
Will be constant throughout the layer and due to the assumed 
constant pressure throughout the field, the density and viscosity 
of gas will be also constant throuhgout the field. Then the 
velocity profile will be again that calculated by Blasius for 
incompressible fluid. By this reasoning, the actual velocity 
profile for lerge Mach's mumber when radiation cannot be 
neglected is something between the Blasius profile and that 
shown in Pig. 3. 


The change in the constant l, MA is 


appreciable, but not great. It decreases from 1.328 for 

M =) to 0.975 for M = IO , or about 30 
percent. However, for QL PE 3 the change of the 
constant is very smell, 

Fig. 1.2 also shows that Eq. (1.88) which was 
obtained by using the linear approximation is fairly accurate 
for very high Mach's numbers. 

As examples, consider first a projectile and 


second, a wingless sounding rocket. Taking the diameter of 


-1l;- 


the projectile to be 6 in., the length 2l, in., the velocity 
1500 ft./sec, and the altitude 32,800 ft. (10 km.), then the 
Reynold's Number based on the total length is 7.86 X 10° 
and the speed ratio is 1.52. From Fig 1.2 the skin friction 


coefficient is 


G = (1.286 x 107)/ ,/7.86=0.000159 


Changing the skin friction coefficient (based on the skin 
area) to the drag coefficient (based on the maximum cross- 
| section), on: obtains 

C Da F 0.0055 
The drag coefficient due to wave formation taken from Kent's 
experiments (Ref. 9) is Cy, — 0,190 
Therefore the ratio of skin friction to wave resistance is 
0.0055/0.190 = 0.029. 

However, the ratio is greatly changed in the case 
of the rocket. Taking the dismeter of the rocket to be 9 ina, 
the length 8 ft., and the altitude of flight 50 km. (see 
Appendix) (164,000 ft.), the velocity 3400 ft./sec., then the 
Reynold's Number based on a density ratio at that altitude of 
0.00067 and temperature 25°C. (deduced from data on meteors) 
is 6.14 x 10º, and the speed ratio is 3.00. From Fig. 1,2, 


the skin friction coefficient is 


l = (1.213 x 10 A/a. = 0400360 


aS 


Then 
Cy = 0.125 
2; 


The drag coefficient due to wave formation from Kent's experiment 


(Ref. 9) is 


Therefore, the ratio of skin friction and wave 
resistance is now 0.123/0.100 = 1.23. If the boundary layer 
ia partly turbulent, the ratio will be oven greater. This 
shows clearly the importance of skin friction in the case of a 
slender body moving with high speed in extremely rerified air. 
It also disproves the belief that wave resistance would always 
he the predominating part in the total drag of a body moving 
with a velocity higher than that of sound. The reason under- 
lying this fact can be easily understood when one recalls that 
the wave resistance of a body is approximately directly pro- 
portional to the velocity, while the skin friction is proportional 
to the velocity raised to a power between 1.5 and 2, Therefore, 
the ratio of skin friction to wave resistance increases with the 
speed. With very high velocities end high kinematic viscosity, 
the wave resistance may even be a negligible portion of the 
total drag of the body. 

Section (11) 

In order to point out the thermodynamic aspect of 
the problem two cases will be considered; the flow of a hot 
fluid along a surface which is kept at a constant temperature 


inferior to that of the fluid, and the case of a hot wall 


-16- 


cooled by a fluid of lower temperature. The problems treated 

in this part have been discussed before in two very interesting 
papers by L. Crocco (Ref. 6). He especially gives an elegant 
treatment of the cooling problem in the case of very high 
velocities ("Hyperaviation"). The author feels that his treat- 
ment is somewhat more general and extended than Crocco's previous 
analysis. 

An interesting general relation between the heat 
transferred through the wall and the frictional drag can be 
obtained using the assumption that Prandti's number, i.e., the 
ratio A/A » is equal to unity. The same 
assumption was used also in the previous calculations. It is 
remarkable that the relation holds also as well for laminar as 
for turbulent flow. The heat flow 4 per unit time and unit 
area of the wall surface is 

= hij ( 97/24 ur 


and the frictional drag z per unit area is 
= EVH 
y = My l Ply 24), 


Using Eq. (1.4) the ratio q /T can be calculated from 


the relation 


E 
l du LEY Ea) 


g Myr ( l. 14) 


where EA is the absolute temperature, and U the 


velocity of the fluid in the free stream, 7), the absolute 
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temperature at the well, Ax and Aur are the heat con- 
duction and viscosity coefficients of the fluid corresponding 
to the wall temperature and Ld denotes Mach's number, 

- Substituting df eo one obtains from Eq. (1.11) 


LL Gwe. Gl) 


= 


ES TT zy aa) 


This is the relation known as Prandtl's or G. I. Taylor's 
formula, first discovered by O. Reynolds. Hence Eq. (1.11) 
gives the correction of this result for compressibility effects. 
In the case A > Fi i.e., when the wall 
is colder than the free etreem, the offeot of compressibility 
is to increase the heat transferred through the wall, However, 
it would be erroneous to interpret this result as an improve- 
mont" in cooling because at high speed the heat produced in the 
boundary layer is of the same order as the heat transferred 
through the wall, In order to determine the efficiency of the 
covling a complete heat bDalence must be made. For this purpose 
Eq» (1.14) does not give sufficient information and the velocity 
and the temperature distribution in the boundary layer must be 
- computed, Such calculations were carried out for the particular 
assumption ay = T i.e., for the particular case 


in which the absolute temperature of the wall is kept constant 


-18- 


at a value equel to one-fourth of the Paran of the hot 
fluid. With the same assumption for the variation of A as 
in Section (I), the results shown in Fig. 1.2 and Fig. lol, were 
obtained. The variation of G JR with Está is 
similar to that obtained in the case without heat conduction 
through the wall, Aloo tho highoot tomporature in the boundary 
layer is very high for extreme Mach numbers. However, the 
temperature maximum occurs some distance from the wall, 

The heat transferred from the boundery layer to the 
wall can be calculated as follows: 


By means of Eq. (1.128), one has 


Hence es 


Lo UA, E Ja, 


f bia a ete) ja 
Syl WE SRL HELO 


h 
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Therefore, combining the above two equations, 


RC 
ley - Lj Es 


(1.16) 


Using Eq. (le7a) and substituting Eq. (1.16) into Eq. (1.5), 


then 


IT) 2 K AR 
Gry, "LE fat Aaa 


where IG == (4º “o /9 + SUN G , aS given 


in Table 1,2, Therefore, the heat transferred to a strip of 


unit width of the wall of length L per unit time is equal to 
L 
/ dee AOA ee 
= / A = pars 
Do LO), Ae bad fe 
i Ay Ta WA 


Now the increase in heat content of the gas per unit time by 


flowing from A=O to X= [ cante calculated as 


PO) 


Q= f [909 im- Ty) 3867 [BZ for 
Yes 
E dr (TH » LAF IES -rik 
(1.19) 


Al A) GAB St TE = KA, TIT 


The viscous dissipation of gas in the boundary layer of unit 


width plate per unit time is thus 
G = Q T a, (1.20) 


Table 1.2 


OuUlIN Fr O 
WN 
pan 
N 
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The total heat balance at different Mach's numbers 
is shown in Fig. 1.5. The "dissipation" curve represents in 
dimensionless form the heat produced by friction per unit time 
and unit width of the plate. The lower curve shows the increase 
(or decrease) of the heat content per unit time and unit width, 
The difference of the ordinates corresponds to the heat trans- 
ferred through the wall. It is seen that cooling takes place 
for fd ee 2.6. Beyond this limit more heat is produced 
by friction than the amount which can be transferred to the 
wall and, as a matter of fact, the fluid is heated, 

In the case A > EA i.e., when the 
wall is hotter than the free stream, the ratio between the 
heat transfer and the drag decreases with increasing Mach's 
number. This is shown in Fig. 1.6 where the ordinate represents 
the ratio between 4/7 with compressibility effect 
(ascording to Eq. (1.1,)) to 9/7 without somprossibility | 
effect (according) to Eq. (1.15)). The calculation was carried 
out for a gas temperature of - 55º F, and a wall temperature of 
180° F, end 300%, IL is seen thet there is no cooling in the 
former case for 2 > 1,69 and in the latter case for 

Jd > 2.08. However, the decrease of cooling efficiency 
is appreciable even at much lower speeds. This emphasizes the 
benefit of the reduction of the speed of cooling air and the 


relatively poor efficiency of cooling surfaces exposed directly 


PD 


to a high-speed airstream. The curves in Fig. 1.6 being 
derived from Eq. (1.1) apply to laminar as well as to 


turbulent motion. 
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APPENDIX TO PART (1) 


ON THE VALIDITY OF THEORY IN VERY RAREFIED AIR 


The hydrodynsmic equation holds so long as the mean 
free path of the molecules is small in comparison with the 
ias of the boundary layer. For this case the thickness 
of the boundary layer is zero at the nose, however, at a 
distance + of the length of the rocket it already amounts to 
3.2 omo, while the calculated mean free path of the air molecules 
at the altitude considered is about 1.1 x 10” om, Hence it 
appears that even for this case the theory can be safely applied. 
This conclusion is substantiated by the experimental results of 
H. Ebert in "Darstellung der Stromungsvorgange von Gasen bei 
néjári gen Drucken mittels Reynoldsscher Zahlen", Zeitschrift fur 


Physik, Bde 85, Se 561-564, 1933» 
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PART (11) 
SUPERSONIC FLOW OVER AN INCLINED 


BODY OF REVOLUTION 


The aerodynamic forces acting on a projectile can 
be divided into three parts: the resistance or drag in the 
direction of the axis of the body, the lift in the direction 
perpendicular to the axis of the body, and the forces due to 
the rotation of tha hady (Magms affect). The first somponent, 
the resistance, is, of course, the most important one, because 
it is the predominating factor in determining the range of the 
projectile. However, in the case of an actual projectile, 
inclination and rotation are always present, and therefore, 
accurate calculation of range is impossible without considering 
the second and third components of aerodynamic forces, 1.0», 
the lift and the forces due to rotation of the body, It is 
found that the solution of von Karman and Moore (Ref. 1) for 
the linearized hydrodynamical equation of axial flew over a 
slender body of revolution can easily be generalized to the 
cason in which the projectile is inclined to the flight path. 
Strictly speaking, the solution is applicable only to a very 
slender body inclined at a small angle to the flight path, 
because second order quantities of the disturbance due to the 


presence of the body are neglected. However, for the case of 


«bm 


axial flow over a cons, von Karman-Moore's first approximation 
(Ref. 1) differs very little from the exact solution of Taylor 
and Maceoll (Ref. 2) for vortex angles up to 10°. Tharefora, 
it is expected that the first approximation of the lift force 
as obtained in this paper can be applied to a pointed 
projectile with fair uucuracye 

If Ø is the potential of the small disturbance v 
velocity due to the presence of a body of revolution whose 
axis coincides with the x-axis, then the linearized equation 
of motion of compressible fluids in cylindrical co-ordinates 
x, r, and O is 


¿LA a 1 RÉ 
2 m 7 gat 7 2? OA? (1.1) 


(/ pe 2 z 
In this equation, V is the velocity of the undisturbed flow 
for which the velocity of sound is sc. If the direction of 
the undisturbed flow coincides with the axis of the body, 
then Ø is independent of O, and Fa. (1.1) reduced to 


Re LES 7% _ 


ax? aL 


(1.2) 


The sulution of this equgtlon when the velocity 


of the undisturbed flow is greater than the velocity of 
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sound, is the same as that for a two-dimensional wave diverging 
from a center, It was obtained by Levi-Civita (Ref. 3) and by 

He Lamb (Ref. 4). Von Karman and Moore (Ref. 1) applied it to 

the present case and showed that it can be expressed as a 


source distribution given by the potential 


a) 


é = y 7 f(x- An cashu) da (3) 


cosh Pe ET) 


r à 
where g= We )-/ . Analogy with a similar case of 


flow of an incompressible fluid leads one to expect the 
solution of Eq. (1.1) to be a doublet distribution given by 
the potential 

o 


Ê = - acoso / SEL “x cosh U)cosh t da (els) 


-/ 
cash E + 


This can be shown to be true, because, if the solution of 


Eq. (1,1) is of the form 


B= cd Alta) 
then Eq. (4.1) reduces to 


Lge OE PL oh ep 
é a? aa Gar A? (1.10) 


Differentiation of Eq. (1,2) with respect to r gives 
fe y? 12 /2% 
+) IX EE, RN E) AER 


Comparing Eq. (b.la) with Eq. (L,2a), it is éastiy seen that 
Eq» (Mel) is a sulution of Bg. (+l); The function f has 
to be determined by the boundary condition 

r 


= 


; 1) 7 e / FIA Ap co ajeak 2 ÁL (4,5) 
Er tag 192 

where Yo is the normal component of the velocity V of the 

undisturbed flow, and R is the radius of the body. 

The complete solution of flow over an inclined body 
of revolution is then obtained by superimposing a cross-flow 
upon an axial flow, i.e. 

À = 
This solution was also suggested independently by C, Ferrari 
(Ref. 5). 

From the velocity potential g, one can calculate 
the pressure distribution over the body and then the aero- 
dynamic forces. However, since the theory is based upon the 
linearized equation, the cross preduct terms of derivatives 


of J, and fo in the pressure calculation cen be neglected. 
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Therefore the following simplification results: the resistance 
or drag can be calculated from the axial flow alone and the 
lift can be separately oomputed from the cross flow. Since 
the resistance was calculated before, (Ref,1), the following 
treatment is concerned only with the lift force. The lift 
ee in a direction perpendicular to the axis of the body 


end the moment about the vortex are thus 


TE po x o 
$ = / T Ap ndé caddy apy) e ity ab doar . 
200 0 


ou YA o0 6 
MT Ap dd tos IX a 2pY [ea cost dba ae 
ia “a YA 2 4 J aX 


whero ^) f is the difference between the pressure at the 
surface of the body and that of the undisturbed flow and f 
is the density of the fluid in the undisturbed flow. 

Eq. (1.5) is a non-homogeneous linear integral 
equation in f which does not have a general solution of 
simple form. However, it is interesting to see how Eq. (2.5) 
simplified in the limiting case when the radius of body 
approaches zeros It is convenient here to use E = Lancha 


as the independent variable, themEq. (1.5) becomes 


LAR tf E: AXAR 
/ ELE MA ENA, ; 
A = “po rá = ob [Lge PME 
2 


uf 


Y xe? 
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Integrating by parts, one has 
A AR q AR 
Y > ýt tera Ey] + FDA 
a 


0 


Now, if tha projactile is pointed nose, the doublet 
strength mst be zero at the nose =O |, thus Ho) ol) 
Let A -70 , ond writing Y instead of Ea in the integrand 


the above equalivn reduces to 


¥ 
A 7 Bi ‘Led 
y e pá Ait) dx Jas: 
J 


Since the cross-sectional area of the body of revolution is 


S = th? » Eq. (15a) can be written as 


x 
TE 
iZ: E EA ee AX 
0 as 4 
or 
x 

Lr 
“gs g Z 

Differentiating, one arrives at 


Fix) = + A (2.7) 


-31- 


In order to calculete the lift, one has first to 


find the axial component of disturbance velocity. Thus 


ag 
LE) --Kesd f AR cosh cash a de 
Gr) 1 [x-AR cor 


cosh 7 a 


= cosg Se “POMBE 
e % N leg) «e? 


£4 
Cos Y 4, eag eS cae AS 
x a L(g) dg = “ib y= Gott ds 


Substituting into Eq. (1.6), the lift force is obtained as 


T © 
L=2pv.f f OD AS sor = EVGA, 


T AX 


Where A, = area of the base section of the body. 


Hence the lift coefficient can be evaluated as 


A 
el ie 
Zt 
in which yo = angle of attack of the body. 


The moment arm d, 1.8., the distance between the 


point of application of the resultant lift force end the 


vertex can be obtained by dividing the moment computed from 
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Eq. (L.6) by the lift force, and thus 


d= (Af) Z is 


mese Am = area of the mean section of the body, 
i.e., the volume of the body divideá by its length, L . 

The results of Eq. (1,8) and Eq. (1.9) are 
identical to those found in Munk's theory of airships (Ref. 6). 
At first sight, this might be surprising. However, if the 
radius of the body approaches zero as assumed, the cross-flow 
pattern is the seme as that for an infinitely long circular 
cylinder moving with its axes perpendicular to the flow. 
Therefore, in every plane perpendicular to the axis of the 
body, the flow can be considered as two dimensional, i.e., it 
is independent of the variable { e Hence Eq. (L,1) 
reducer simply to 


TAO O g 


In? YA IA (Le lb) 


This is immediately recognized as the equation of motion for 
two dimensional flow of incompressible fluids, which is the 
basis of Munk's theory. 

Due to this two dimensional character of the flow, 
the distribution of doublets is not affected by the change in 


Mach's number, which is only comected with the independent 


¿Je 


variable x , and, therefore, the lift coefficient and the 
moment arm are also independent of Mach's number as shown by 

Eqs (1.8) end Eg. (269). This oan aloo bo soon from tho fact 
that when r approaches zero, the variable Eaf- Xh ahu 
—» X and thus the effect of oœ , which is a function of 
Mach’ s number, is removed. To study the effect of Much's 
number on the lift of the body, one has to go back to Bq. (2.5). 
To avoid the difficulty of solving this integral equation, the 
"indirect method" of solution can be employed, l.0., take a func- 
tion f and determine the necessary shape of the body ta comply 
with this function f. 


Taking the simplest case 


flx- th chu) = K (2-an cosh a) 


where K = a constant, Then 
0 
Ê, = -Ad csp pa (X- AN Comu) Cosh u du 
tosh E 
E 2 AR 
Kacoso) + EJ? = F cosh A 


ånd the boundary condition reduces to 


- hat fã 2 ya 
= E Malal! + th a 
Therefore, the solution, Eq. (2./ja), is evidently a solution 


for a cone with half vertex anglo € , if cot € = yA . 


By putting cof é = E the boundary condition can 
x 


-3) po 


be written in the form 


2 
a E CAE, + cah “ef (2.5b) 


For any given value of vertex angle and Mach's number, the 
corresponding value of K can be obteined from Eq. (Z.5b) 

In order to calculate the lift, one has first to 
find the axial component of disturbance velocity. Thus from 
Eqe (Dela) 


ELA seo Ra 


IX hep =A “E AJ , Waz a 
= wh cos o [627 


Substituting into Eq. (1.6), the lift force is found to be 


= 2E Va, i 2% 
= att EE a o) 


¿TR AX 


= AN fst LV cos € 
AA ie 
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where e = lateral surface area of the cone. Therefore 


Pe ee ties a a o cot EP 


But from Eq. (1.5b), t is obtained as 


B Lu 
a? yea + coah Cf 
Hence y = A ar : (2,10) 


Z 
[S—/ 
where A = Loy ST 
-/ 
SASE + oh S 
In the limiting ease when € approached zero, 
fy AL 
which agrees with Eq. (1.8). Similarly from Eq. (1.6) the 


moment coefficient is 


Ĉ, = moment about vertex — _ 2 
m LÍO ça 


VAL 


which satisfies Eq. (1.9). 
Both Eq. (1.8) and Eq. (2.10) show that the lift at a given 


Mach's number is proportional to the angle of attack of the 
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body. This is a general characteristio of flow around a body 
without separatione If Une fluid separates from the Lody mud 
creates a "dead water" region on lee side of the body, then 
the lift will be proportional to the square of the angle of 
attack as was shown by W. Bollay (Ref. 7). The problem 
whether the fluid separates or not can only be answered hy 
experiments. From the experimental data now available, (Ref. 
8), it seems that the flow is continuous without separation, 
end, therefore, the lift is proportional to the angle of attack 
of the body. 

Fig. %e2 is the result of computation using Eq. 
(1.10). Calculations were carried out for values of A. Zi, 


because the value of “Æ = / corresponds to REA 


/ 3 


where É is the wave angle. For A < 7 it is found 
thet P < € - This means thet the wave angle is maller 
than the vertex angle which is, of course, impossible, There- 
fore, A, = / marks the limit of validity of this solution, 
In fact, even when A, is near to 1, the solution must be 
considered es qualtitave only, since in this region the effect 
of the surface of the body on the shock wave cannot be 
neglected. 

To generalize the solution for a body of revolution 
with a sharp point at the origin and cylindrical shape at 
infinity, it is simplest to use a step-wise doublet distribu- 


tion. Consider the points Fas A Po... Fe de E, of the 


=37- 


meridian line of the body, end designate their co-ordinates by 
A, LA y A, ko, A Ay Ra te A, Lu and the correspond- 


ing values of y- «Æ by A E, E ores F 


2 


Then the boundary condition of Eq. (2.5) can be written as 


~ cos ee 


This condition actually gives a set of N equations to determine 
the N constants A; + This set of equations can be 
solved rather easily because each following equation in the set 
only contains one more which does not appear in the preceding 
one of the set. When KA $ are determined, the lift can be 
calculated by using Eq. (1.6). The pressure over each section 
of the conical surface is constant. The lift and moment co- 


efficients are thus obtained as 


yp = = 
> E o Cn a eat) Fie Fay [En Bins)" MES “ff 


(a 
7 tato ori the Oy a fee 24 ph ha He ) 


EE -7 AE ‘If (2.12 


-38= 


where E is the last point on the meridian line, Ha is 
the base radius and / is length of the body. 

Fig. 2.3 is the result of the calculation using Eq. 
(2,12) for a body of revolution with "6-cliber head" and 
total length of 4.8 calibers, when it is travelling with a 
velocity 2,60 timos that of sound ( a = 2,5). The lift 
coefficient is considerably higher than that of a cone at the 
same angle of attack and at the same Mach's number, evidently 
dus to Lhe cylindrical part of the body, The position of the 
resulting lift force is also shown in the figure. Since, as 
mentioned before, lift and drag are independent in first order 
approximation, the calculated lift coefficient can be combined 
with the drag coefficient taken from experiments and thus give 
some information on the magnitude and direction of the result- 
ing force. Fig. 2.) shows the method applied to this projectile 
with the drag coefficient taken from Kent's experiment (Ref. 9). 

If the projectile has a length of l.Z3l diameters in- 
stead of 1.8 diameters and has the same nose shape, and, if its 
center of gravity is located at a point 2,68 diameters back of 
the nose, then the calculated moment about the center of 
gravity can be expressed as lids = S VA, “yo AL 
where FL) = 735 for + = 269 . 
This compares closely with the value AL) = /OF 


extrapoleted from R. H, Fowler's experiment (Ref. 8), for a 
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projectile of the same proportions, This shows that the theory 
developed in this paper can be applied to a projectile with 


fair accuracy. 
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APPLICATION OF TECHAPLIGIN'S TRANSFORMATION TO 


TWO DIMENSIONAL SUBSONIC FLOW 


“The equations of two dimensional irrotational 
motion of compressible fluids, assuming that the pressure 
is a single-valved function of density only, can be reduced 
to a single non-linear equation of the velocity potential. 
In the smpersonie case, the problem is solved by Prandtl, 
Meyer and Busemann by means of the powerful method of 
cheraoteristios. The essential difficulty of this problem 
lios in tho subsonic case especially when the velocity 15 
near to the velocity of sound. The first logical step is to 
linearize the equation based on the argument that the disturb- 
ance super-imposed on the parallel rectilinear flow due to the 
presence of a solid body is sufficiently small compared with 
parallel flow, This makes the second and higher order terms 
of disturbance potential to be negligible. An exemple of 
this method is the well-known theory of thin airfoil due to 
Prandtl and Glauert. But the presence of stagnation point at 
the nose of the airfoil makes the application of the linearized 
theory questionable, at least near this region; because thara 
the disturbance due to the presence of the body is no longer 


small, On the same ground, the theory breaks down in case of 


h2- 


bodies whose dimension across the stream is not small compared 
with the dimension parallel to the stream. The next method is 
thet derived originally by Janzen and Lord Rayleigh. They 
solved the equation by successive approximations. However, 

the process is very tedious and the method convergent very slowly 
if the velocity approaches that of sound. 

Molonbroock (Ref. 1) and Tochapligin (Ref. 2) 
suggested the use of the magnitude of velocity and inclination 
of velocity to the x-axis as independent variables and were 
able thus to reduce the equation of velocity potential to a 
linear equation. This equation was solved by Tschapligin (Ref. 
2) and recently put in a more convenient form by F. Clsuser end 
M. Clauser (Rof. 3). The solution is essentially a series 
each term of which is a product of hypergeometric function and 
circular function. The chief difficulty in practical application 
of this solution is to obtain a proper set of boundary conditions 
in the transformed plane, or the hodograph plans. 

Tschapligin (Ref. 2) showed that a great simplifica- 
tion of the equation in hodograph plane results if the ratio of 
specific heats of the gas is equal to -1. Then the equation 
becomes the equation of minimal surface whose solution is well- 
Inova, However, at first, the hypothetical value of ratio of 
specific heats (all real gas has the value for this ratio 


ranging from 1,00 to 2,00) makes the practical application of 


“igs 


Tschapligin's theory questionable. It was Demtschenko (Ref. 1) 
and Busemann (Ref. 5) who made the meaning of this special velue 
of ratio clear. They found that this special value of ratio of 
specific heats really corresponds to take the tangent of -v 
curve of gas instead $- vv curve itself. However, they 
limit themselves to use the tangent at the state of rest of the 
gas. Thus their theory can only apply to velocities up to 0.5 
times sound velootiy. In this part , the theory is generalized 
to use the tangent at the state of gas corresponding to the 
undisturbed parallel flow. Therefore the range of usefulness 
of the theory is greatly extended. In the first section, the 
general theory will be developed. In the second section, the 
theory will be applied to the case of symmetrical Joukowsky 
airfoil at zero angle of attack. 

Section (I) 

If $ is the pressure and 147 is the density of 
gas, the adiabatic process is expressed as a curve in the pa 
plane as shown in Fig. 3.1. Now conditions near to the point 

bi, v, can be approximated by the tangent at this 
point, The equation of the tangent at this point oan be 


written as 


P,- = l- v)= elg 9”) (3-1) 


where Ç is the density of the gas. Now the slope É 


= 


must be equal to the slope of the curve at the point PY 
- sdb) _ ¡de de Aliada Ei 
C= Gy) = Cde dv), = - (FE) se as 


2 
Therefore C E a ¢ 
Hence the approximate p- $ relation near $, 7 5, 


can be written as 
22 / f / 
ds AS, de 7 (3.2) 
From the generalized Bernoulli's theorem, the 
following relation is obtained 


? 
Lopra See oe ap 
7 Wa ss ai pe $ 


= 


(3.3) 


where 4/7 is the velocity of gas, the subscripts 2 3 
indicate two conditions. But from Eqe (3.2), ? can be 


expressed as a function of $ » thus 


Ap = 4, Ss de (3.0 


Substituting into the integrand in Eq. (3.3) and integrate, 


the following relation is obtained 


DA eee A ie 
fa- fats [ge pf 


Now if wy =o, We = ur DEY and 5 =ç + then 


4 
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>` 
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ge e (3.5) 


where the subscript O denotes the rest state of the gas. 
If the sound velocity 4 is defined as the derivative of 


$ with respect to E- , then Eq. (3.4) gives 


sê e ati 95% constant (3,6) 


Therefore Eq. (3.5) can be written as 


$, ar? 
(E) =/- 35 


or 
g 2 
7 (- E (3.7) 
Furthermore, from Eq» (3.6), Sa” = Sa , 


thus Eq. (3.7) can also be written as 


e ja LL (3.8) 


It ia interasting at this stage, to notica thet 
from Eq. (3.8), the density decreases as velocity we 
increases, as it is expected. Thus from Eq. (3.7), the 


velocity of sound of the gas will increase as the velucity is 


-y6- 


increased. This is just opposite to real gas, because in the 
case of an adiabatic flow, it is well-known that the temperature 
of gas decreases as the velocity of gas is increased and thus 
the sound velocity also decreases. However from Eq. (3.7), the 
ratio A or Mach's number increases as the velocity um 
increases. But this ratio only reaches the value unity when 

ur = œ , or from Eq. (3.8) when $ =0J + Itis 
thus seen that the entire region of flow is subsonic and thus 
the equation of motion is always of elliptic type. This may be 
considered as the physical reason why the complex representation 
of velocity potential and stream function is possible in all 
cases, as will be shown in following pages. However one should 
bear in mind that the portion of tangent that could be used as 
an approximation to the true adiabatic curve, must lie within 
the first quadrant. Thus the upper limit of velocity for 


practical application of the theory is when a =P « By 


using Eqs. (3.2), (3.7) and (3.8), this upper limit is found to 


be 
MEN en 
( Ar) a 
"max 
Or by putting 4 = pee , the above equation 


reduces to 
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tar 


The values of (yyy for different values of SÉ 
7 á / 


are shown in Table 3.1. 


Table 3.1 
= ar az (a, ree 
0 o9 2.186 
0.2 10.91 2.195 
0.4 5.56 2.205 
0.6 3.78 2.265 
0,8 2.92 2.335 
1.0 2.405 Belos 


It is thus seen thet for most applications of this 
theory, A will remain positive. However due to large 


deviation from the true adiabatic process at high values of 


uw AN 
w , one has probably to limit the ratio Cos to 
about 2. 
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Now if the flow is irrotational, there exists a 
velocity potential $ such that 


bey > 


where 4, Y are the J Y components of the 
velocity 4 . To satisfy the equation of continuity, the 


stream function p is introduced. It is defined by 


£ya i’ Euro BE 
RR g OS Sy (3.10) 


Now if the angle of inclination of the velocity 47” to the 


x-axis is Ê , then from Egs. (3.9) end (3.10), one has 


Ag = y cos p dx + ur sing dy 


Ap = -u é MA d+ a & cos p dy (3.11) 
Solving for XY and dy , 
- inf Ad 7 Sng Edy 
(3.12) 


— æn as E 
dy = Pe db + ms LV 


So long as the correspondence between the physical plane and 
HX, Y) 
the hodopl i t th ticall: 2 
o prono Ss ono o ONG, or me come ca y au, tr) A 
one can express x P y as functions of 44 Z and so 


also 4 and YP as functions ‘4 and f e Thus 
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dé = di dut by Ap 
A= Vo dur + “as Ap (5.15) 


where primes indicate the derivative, and subscript indicate 
the variables with respect to which the functions are 
differentiated. Now substitute Eq. (3.13) into Eq. (3.12), 


the following relations are obtained 
/ cos p g! Mob Saro wu CAES Linb Bye 
dam [Pty PPV lO LEV Me 


(3.14) 

[Ep 4 cos Sb Ss, OB & apf 
tya E het EE Ve der Gat Feu 
Since the left-hand sidesof Eq. (3.14) are exact differential, 


one can apply the reciprocity relation and obtains 


ip (SEG, EE ee ay May) 
Po. 15) 
+ Sin / 1 . 
(EL, A COS Y) - E (A, oh ey) 
Carrying out the differentiation, and cancelling identical 


terms in left-hand and right-hand side, 


— wn / _ OSB Sain’ - CB F’, SB Sy Misao! 

ee Pe eG ee A 
3.1 

E aa fm LA f 2 / Cos A Na 5 
u AEV =- Deg- E 
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Using Eq. (3.7), Eq. (3.16) becomes 
“sn / e EA Z LOA NB S sip 
Ww a IA 5 Ya = w’ mid é “ So Y 


COS “So / y 4 + (3617) 
y tab W =~ Ed, Sb 


wr TÉR A wh 


/ A / lo? 
As in both equation. g 4 and VA 
quata ar j q E , A 
are connected with a proportional factor, one can solve for 


them, and 


o es ; 
wo a ie 


> 


(3 +18) 


Now Eq. (3.18) can be further reduced if a new 


variable (i) is introduced. (J is defined as 


E, w (3.19) 


a 
i 
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/ = 7 (3.20) 
ed 


This is the fundamental set of equations for the present theory. 
It can be easily recognized as the Cauchy-Riemann differential 
equation, and thus g + iP must be an analytic 
function of A +: . However for the convenience 
of numerical calculation, a new variable WY is used instead 


of (0  , such that 


W = 4, . e (3.218) 
Or by integrating Eq. (3.19), 
24,4 


W = 3.21 
/ y atu” + 4, q 


Hence by inverting, 


4 Ay W (3.22) 


Uy = 
ha? - w? 


Thus by substituting into Ege (308), 


Do haf + Wwe — 
§ Z 2 ê 
La, - W 


if another set of new variables U - w cosp 


and V = W sm B are used as independent variables, one 
has 


JU + Ii = W jaspè a 
do = Wf- yg + cap f 


Using Eq. (3.24), Eq. (3.20) can be written as 


espe a : 
esp dE 4 sm 626 = uma BE - cos ze 


Ú 


— Sin À ZA ae Y 
sing 26 7 tos BOG op y sng 2 


These equations aro satisfied by 


_ 26 ad . ap 
Y AV), a ACV "QU (3.25) 


These are the Caushy-Riemann differential equations, therefore 
the complex potential / = p+ ¿yo is a function 
of O - ¿y = W | e Or. 

Pr iP = EUA) = AW) 
Hence A e yp E F (U+iy)- Flw) (3,26) 
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To transform from hodograph plane back to physical 
plane, the expression of Y end Y in terms of ZA and 
V must be found. By using Eqs. (3.22) end (3.23), Eq. (3.12) 
can be written as 


Var we 
A pe] E aro ha? 


y= fr; je te UM fre fof 


2 2 2 
where W = +y . These equations can be combined 


into one equation by means of Eq. (3.26). Thus 


de = duvidy= dE web 


“Rai (3.27) 


For practical application of the theory to tha fl ow 
over an obstacle, the computation proceeds as follows: (1) 
Find the complex potential for the flow of incompressible 


fluid over the obstacle, sey 
w, GCE tig) - w, GE) 


where W, , is the velocity of parallel rectilinear undisturbed 
flow, end 3 , n the space coordinate of the physical plane. 


(2) nowt =F = W GIS) . Here jf isthe 


suis 


transformed undisturbed velocity, to be interpreted as Eq. 
(3.21). But the complex variable É has no direct 
physical meaning. (3) Using the above value of F , Bae 


(3.27) can be written as 
LIME IAE Y? po 
de= dí - A E) (iz) 
Integrating, 


Z= $- LES ÉS de (3.25) 


Thus it is seen that the complex coordinate in the physical 
plane of compressible fluid is equal to the corresponding 
complex coordinate in the physical plane of incompressible 
fluid plus a correction term. The factor before the integral 
depends upon the Mach's number of the undisturbed flow only. 


By using Eqs. (3.7) and (3.8), and (3.21), one has 


(Lf 


aoe Pra Jr Eu e 


The integration constant of the integral in Eq. 


(3.28) is not important, because it only means a translation 
of the whole 2 - plano. (hj Tho velocity w 


corresponds to 2 can be computed by starting with 


T AAE e ie 
We fe Maen UV 


By means of Eq. (3.29) 
ee 


Wee LOLA 
nfl! FOG) 


Thus by putting AO — tu, , one obtains 


M E 
Thus Lyr (3.30) 
LAY TS ET 


Using Eqs. (3.30 and (3.29), the ratio a can be 


~~ 


calculated easily. (5) To determine the pressure acting on 
the surface of the body one has to use Eq (3.2). With some 


manipuletion, the following relation is obtained: 


2 
43 m7 > (5) ae FE) 


as) 
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But 
e) 


(3.31) 


Therefore 
(3.32) 
f Pp; 2 / tr? “yf 
= tay IAEA EG 
Fsm EF NG NA 


Section (II) 


In this section, the general theory developed in 
Section (I) is applied to the simple case of flow over a sym- 
metrical Joukowsky airfoil at zero angle of attack. The 


complex potential in the circle-plane (see Fig, 3.2) is known 


to be 
-b)+ 4 
u fer) + | 
i ¡A Ue (3.33) 
where @& = radius of the airfoil circle, ge 


eccentricity of the airfoil circle. The relation between the 
airfoil plane and circle-pjlane is the well-known Joukowsky 


transformation it 
OST 
(3.31) 
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if the radius of the transforming circle is unity. 
Now the starting point of the calculation is the 


function to find MW AE 


= 2 
dE dae (E) 
WdF = Fe di ay Mp 


Therefore yf = n ji- Ui pf H ei Ee AN 


Thus the correction term in Eq. (3.28) is 


ig JWA = AU rd] (3.35) 
where Ios [fr a q) 
dl 2 
2 [lr SD 


PE E tee de 
5 AL [lps E arte 


These integrals can be easily computed and simplified, noting 
thet a-5= / . If (ub )=ne4 and 


/ =f -268 “a 
Ll, = 2] A - CS el a ge Laplace] 
LF i t-A log al eh E) dla A lag ae 
| 2) 00 qe 28 -38 .36 
AA Ur A)e M+ É e fe f (3.36) 
Separating the real and imaginary parts, and adding, 
el (1 O a FA od - + cos 38) 
2 
F A / LE hoy [nf + Es AFA paar 


acs (1-A*) cos YY 4 F 12A) cm 38 f 
(3.37) 


Im (144-4) = a(o s08 +$ SIE) 
+ El f (424-23) sing - ff 1-22) sm 2d 


“+ lita) sab (AT by EE Haaf 


vA 
EE 5? 


These give the correction term to X and Y 


coordinates e 


The transformed velocity MW over the surface of 
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the airfoil can be easily found by means of graphical method 
(Ref. 6). Then the actual velocity and pressure can be computed 
by using Eqs. (3.22) and (3.32). 
Fig. 3.3 shows the result of caloulation for the case 

4 = 1.20 and om = 0.550, The nose of the 
airfoil is somewhat rounded by transforming into the case of 
compressible fluid. However, the pressure gradient is steeper, 
as would be expected. The main defect of this type of cal- 
culation is that during the transformation from incompressible 
flow to compressible flow, the shape of the body is also 
changed. To isolate the effect of compressibility of the 
fluid, it is necessary to bring back the original shape of the 
body. This is done by first deforming the original Joukowsky 
airfoil, such that the final profile after correction for 
compressibility is same as the original Joukowsky airfoil. 
The amount of deformation is obtained from the calculation 
assuming that the airfoil was a Joukowsky airfoil at start. 
That is the effect of deformation on the correction term of 
Eqe (3.28) io nogleotod. This is allowable because the 
quantity neglected is a second order quantity. 

This deformation of Joukowsky airfoil can be 

carried out by using the method developed ty vou Karman and 
Trefftz (Ref. 7). However, for some practical reasons, the 


Karman-Trefftz method is somewhat modified; 


-bla 


Fig. 34a shows two eirfoils, o having the same 
chord, one is Joukowsky airfoil desired and the other is the 
airfoil resulted from the first step calculations Now apply 
the Joukowsky transformation to this figure, then the Joukowsky 
airfoil will become a circle 0 while the other airfoil 
a near-circular shape, asshown in Fig. 3.lb. The desired 
deformed Joukowsky airfoil will appear like C 2 in this 
figure, The difference between É, and C is just 


opposite and equal to that between E , and C « Now let 


C; be written as G=% e? e Obviously, 


l= St KA (3.38) 


where ao) wil be small compared with 1. The function 
which established the conformal transformation of the outside 


of this boundary Cc into outside of the circle É / may be 


denoted by 

Z = 

5, = Sa LIFA (3.39) 
where C E E 2 have their origin at the center 
and the absolute value of Ll E ) is again small 


compared with l. Then it is shown (Ref. 7) that 
A 2) + fe A H Ss) ee (3.10) 


In order to calculat e + ‘a 2 ) wo develop the 
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function gl 2 ) in a Fourier series: 


oo 
(8) = Gy cos NG 
HO Zo Gula) 
Here only consine terms appear because the airfoil is 


symmetrical about the chord. One the other hand, the complex 


function PEL) has the form, for [$2 [> f : 
> e 
4 A 
+ é Cr) = E? (3.42) 
ô 2 
Now put Es =x ¿4 , then (3.40) is satisfied by 
| Cy = -åy 
Thus 


Mess a (3.43) 
a G 


It can be easily seen that the velocity around the deformed 


Joukowsky airfoil can be calculated as 


AC (ehl) 
where Wr = velocity around the Joukowsky airfoil. 


Now trom Eqs. (5.59 amd (3.43), 
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= {14 Žiada Ad - i Dil) agin ad 


Neglecting small qum tities of second order, end noting Eq. 
(3-42), 


|} = /+ Zi (2-1) 4, cos nb 


Do co 
= ee a Nay cos NO -Flyin n (3.45) 
0 O 


oo 


=/4 Zn sin nf - 28) 


A trial calculation shows that the convergence of 
the coefficients ay is not very good. Therefore, one 
must avoid manipulation on the Fourier series as required by 

eo 
Eq. (3.45). This is possible because — 7 Ay Sin nO 
is known to mathematicians as the allied or conjugate series 


oo 
of a ay Cos nd = - It io also known 
O 
(Ref. 8) that if pb) = a Gy COS VB 
ô 


T 


a TREM HOP AI) 
a = f E df 


Therefore 


YA 


i tE) -lA 
D G Snn = a | ORR y 


Integrating by parts, 


“+ SUEDIA u jo y [terns dei 


a El "y, Hars)~ p+ SJE) 0 
aaa 


Hence Eq. (3.45) can be written as 


[tile |g op [La pop ge tH, 

Fiz)" /- cos E lo) +146) | 
The integral is evidently convergent for any continuous regular 
function 4/0 ) , because then the integrand is always 
finite. Its evaluation can be done numerically. 

Figo 505 is the result of calculation for a 

Joukowsky airfoil with the thickness paremeter b =0 20 g 
at two speeds, M = 0.450 and 0.550. The suction 


a, 
peaks are considerably higher with higher speeds. Also the 
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positions of pressure peaks tend to move backward with increas- 
ing speed. Both are in agreement with the experimental results 


obtained by J. Stack ( Ref. 9). The values of ( PAYA ga? 


“for M = 0.550 end SL = 0,450 at which 
f 4 


real air will attain a velocity equal to the local sound velocity 
are equal to -1.653 and -2.755 respectively. It is thus seen 
that the effect of compressibility on pressure distribution is 
appreciable, even when nowhere the local sound velocity is 
reached. One should, however, bear in mind that the effect on 
the force coefficient of the airfoil will probably not be so 
marked as with the pressure distribution, because the resultant 
forse on the airfoil is the algebraic difference of pressure 


force acting on two sides of the section. 


APPENDIX TO PART (111) 


COMPARISON WITH OTHER METHODS 


Tn order to check the accuracy of the method 
developed in PART (111), the flow over a finite circular 
cylinder with its axis perpendicular to the direction of 
undisturbed flow is studied. The method exposed in Section 
(11) of PART (III) for correction of shape of body is used. 
The following is the result of calculation for velocity at the 
top of the cirouler section, compared with results by other 


methods [collected by E. Pistolesi (Ref. 10)/. 


E = 0.100 
Method 7 at top of section 
Part (111) 2.268 
_ Rayleigh 2.206 
Poggi 2.194 
Taylor's 
Electric 2.188 
Analogy 
A 


w = 2.000 for incompressible fluid, 

t 

Thus the present method gives e higher value. However, the 
flow over a cylinder is rather an extremo case. Because the 


difference between the velocities to be calculated and the 
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and the undisturbed velocity is large, and thus this 


epproximete method involves lerger then usual error. 
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PART (IV) 
FLIGHT ANALYSIS OF A SOUNDING ROCKET WITH SPECIAL 


REFERENCE TO PROPULSION BY SUCCESSIVE IMPULSES. 


Introdution 

In 1919 R. E. Goddard (Ref. 1) published the 
historically importent paper which suggested the use of nitro- 
cellulose powder as a propellant for raising a sounding 
rocket to altitudes beyond the range of sounding balloons. To 
delermine Lhe feasibility of this propellant, a series of exe 
periments had been carried out and it was found that thermal 
efficiency of 50% could be expected if the powder was exploded 
in a properly designed chamber and the resulting gases were 
allowed to escape at high velocity through an expanding nozzle. 
In 1931 R. Tilling used a mixture of potassium chlorate and 
nephthalene as propellant and actually reached an altitude of 
6,600 feet. More recently, L. Damblenc (Re. 2) made static 
tests with a slow burning black powder and from these 
estimated that a height of 10,000 feet could be reached using 
a two-step arrangement. The remita so far reported offer an 
incentive to further enalysis. 

The effect of decreasing gravitational acceleration 
on the maximum height reached by a rocket has been considered 


by A. Bartocei ( Ref. 5). However, he assumes that the 
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rocket itself has a constant acceleration during powered 
flight. L. Breguet and R. Devillers ( Ref. 4) also con- 
sidered the effect of the variation of g, To simplify the 
analysis, they assumed that the acceleration of the rocket was 
equal to a constant multiple of ge Since the sounding rocket 
for practical reasons will be propelled by a nearly constent 
thrust or a uniform rate of successive impulses, in Section 
(II) the euthor has studied the problem anew according to 
this mode of propulsion. 

When the sounding rocket is ascending through the 
eir the maximum height reached is less than that reached for 
flight in vacuo. Recently, studies have been made of the 
problem by W. Ley and H. Schaefer (Ref. 5) and by F. J. 
Malina and A. M, O. Smith (Ref, 6). On the basis of the 
latter study a group of new performance paremeters have been 
isolated from the general performance equation, end these are 
discussed in Section (III). 

Notation 

Referring to Fig. (4.1), the following notation has 

been used throughout the paper: 


Ws weight of propellant and propellant container ejected 
per impulse, lbs, 


k = ratio of container weight to sum of container and 
propellant weight ejected per impulse, 


À= (1-£) 


W wz initial weight of rocket, lbs. 
0 
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M, e initial mass of the rocket, slugs. 


Wp = 
E 


instentaneous weight of rocket, lbs. 


ratio of initial weight of propellants to initial total 
weight of a rocket propelled by coustaul Lhrustl. 


ratio of initial weight of propellents to initial total 
weight of a rocket propelled by successive impulses. 


T AA 


number of impulses per second. 

total number of impulses cccuring during powered flight, 
time, sec. 

interval between impulses, sece 

initial acceleration imparted to rocket, tt/ sec. 


acceleration of gravity at the starting point of flight, 
ft/ sec“, 


acceleration of gravity above the starting point of 
flight, ft/sec“. . 


effective exhaust velocity of ejected propellant, ft/sec. 
instantaneous velocity, ft/sec. 


velocity imparted to rocket by the L +) impulse, 
ft/sec. 


velocity at the end of the nh interval, ft/sec. 


velocity of sound corresponding to the atmospheric con- 
ditions et the starting point of the flight, ft/sec. 


velocity of sound corresponding to the atmospheric con- 
ditions at the height reached by the rocket at the time 
t, £t/ sec, 

Machts number = rly 


veloaity of rocket at start of cossting flight, ft/sec. 
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35 


i 


velocity of rocket at start of coasting flight if ¢ is 
constant and equal to g,» ft/sec. 


altitude above sea level, feet. 


height reached at the beginning of the 224 
interval, feet. 


height reached at the end of the ath interval, feet, 
height traveled during powered flight, feet, 


height travelled during powered flight, if g is constant 
and equel to g,, feet. 


height travelled during coasting flight, feet, 


height travelled during coasting flight, if g is 
constant end equal to go, feet, 


height travelled during powered flight and coasting 
flight, feet, 


= height travelled during powered flight and coasting 


flight, if g is constant and equal to g, feet. 

radius of earth, 2.089 x 108, feet. 

drag on rocket due to air resistance, lbs, 

drag coefficient of rocket shell. 

drag coefficient of rocket shell at the velocity of sound. 
drag-weight factor (discussed in the section on the effect 
of air resistance). 


mass density of air at the starting point of the flight, 
slugs per cu. ft. 


ratio of air densities at altitude and at the starting 
point of the flight. 


absolute temperature of the atmosphere at the height 
reached by the rocket at the time T, °F. 


absolute temperature of atmosphere at the starting point 
of flight, °F. 


“Tee 


A = largest cross-sectionel area of rocket shell, sq.ft. 
d = largest diameter of rocket shell, ft. 
2 = length of rocket shell, ft. 

Section (1) 

An approximate method of calculating the maximum 
height reached by a rocket propelled by powier was developed 
by Re He Goddard (Ref. 1). To simplify the analysis a contin- 
uous loss of mass was assumed and the problem wes so stated 
thet e minimum mass of propellant necessary to lift one pound 
of mass at the end of the flight to any desired height was 
determined. However, if high-powered powder is used, the rate 
of burning is so rapid that the propulsive action is insfanten- 
eous, The rocket is thus acted upon by an impulse rather than 
by a constent thrust, 

In the following analysis, it has, therefore, been 
assumed that the propulsive force is an impulsive force, a 
the force acts for such a brief intervel of time that the 
rocket does not change its position during the application of 
the force, although its velocity and its momentum receive a 
finite change, If the combustion process of the propulsive 
unit takes place at sonstant volume this assumption is 
justified. Further, a study of interior ballistics of small 
arms reveals that the period between the ignition of the 


powder charge and the bullet's arrival at the end of a two- 


«Tha 


foot barrel is of the order of ll; ten-thousandths of a seconds 
If the gases are not restrained and their travel through the 
burning chember end the nozzle is of much shorter length, as 
is the case for the rocket motor, even shorter periods of 
duration of action can be expected. 

Assuming that the propulsive force acts as an im- 
pulse, then the motion of the rocket can be calculated by 
Newton's third law, which states that impulses between two 
bodies arc oqual and opposite. Hence, eating tho momentum 
of the exhaust gases to the momentum imparted to the rocket, 
using the quantities defined in the list of notation and 
referring to Fig. 4.1, the following relation can be written 


for flight in weno 


Aw np = Ma 
e (an 


where A= /-k ana W, = W -Aw (4.2) 


Or Av ES WAC = zte f) (1.3) 


During the intervel between impulses, AL j 
the velocity is reduced by the action of gravity so that at 


the end of the 2% interval, the velocity of the 


-The 


rocket will be 


/ 
Y= Y -qat =U, + A - gat (hek) 
Therefore 
S=. 
v, ¿e AV -AA At 
E> S 4, (1.5) 
JS=/ 
Substituting for 4 Us from Eq. (hel) 
f= 
/ É ‘he 


—— -A9 Al 
Po NG SE 4 (1.6) 


Or wv’ = SAC Saya 
i 


A N o 

s- 
where = 

S=/ E) 


The height gained during each interval will be 
represented by the area under the velocity curve in the 


interval, or 


tha = Yat 12414] (4.7) 


Ga 


Therefore, at the end of the Nth interval which is 


the end of the powered flight, the height will be 


Aen m 5 
Hp, ee Y NÉ + L- g, (at) 
l=1 
Substituting for ae ite value in Eq. (4.6) 
REN z, SER 
PA ea — -agat f+ giat 
A= N Sa Es) 
a =N 
FR 
ERA 
a= ae) 
= GAC Ac a 
ge RE 7 (st)? (1.8) 
AN 
where a = NA 1-2 
Z pe pie i 
nr 


The maximum height reached will be the sum of the 
height at the end of powered flight and the height travelled 


during coasting or 
Fray = H, +H, = Ho E 22% (4.9) 


To calculate the maximum height one has first to 


evaluate the sums S, and Soe Noting that 


my fons 


/ a 7. pi Dm 


sl) do 


S] cen be written in the form 


y =/ 
E a 
Sar (tar 
= -X — 
ye x ax 
g Peço Ms 


Putting X= qe the above integral becomes 
N. A — NY / = € re 
eee E a y 
a E YA € Y e! e? dy 


2 HE) WE a 


(4.10) 


where Y(a)- É [taraf , the so-oelled psi- 


function (Ref. 7,8). 


Similarly, So can be summed as 


= 7 / N-f E nave cz) VE wjf (4.11) 
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Substituting Eqs. (4.10 end (1.11) into Eqs. (4.6) and (4.8), 
and then into Eq. (4.9) finally 


Anas, = ty- e del (E ) Y. Nf (4.12) 


where y = 4 


Zo m p WA) WE) 


For convenience of calculation in Fig. 4.2 the 

quantity Y 
£ 

of 5, . 


is plotted against N for different values 


It cen easily be shown that when N 


W = 


so that Eq. (4.12) reduces to 


1 


2,2 E á Z 
H = (7) | 
MAX, 24 E 


( L . 12a) 


Also, as N —> eo. e = bogli-E/) thus Eq. (4.12) 


reduces to 


oe Es) 
His, = dE [4 dez E ae 7 (4.120) 


+L 


f 


= -LnÃt-G= aa L.12 
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The quantity a | 


a can be considered as the initial 


acceleration of the rocket if N — 20 + Itis 
interesting to notice that Eq. (L.12b) is the equation ob- 
tained by Malina and Smith (Ref. 6) for calculating the 
meximm height of a constant thrust rocket, as expected. 

Fig. 4.3 shows the variation of pine 
with p1c/%, for different valves of €,” 
and for four valves of N. These curves show that when the 
total number of impulses, N, becomes larger then 100, the 
maximum height reached is imperceptibly chenged by increasing 
the number. 

At this point it is necessary to discuss the 
similerity existing between a rocket propelled by successive 
impulses and s rocket propelled by constant thrust. The 
former loses not only the mass of the propellant, but also the 
containers for the individual charges, The difference in | 
effect on the rocket between the propellent and its containers 
is that the propellant has en effective exhaust velocity, 0, 
while the ejected containers leave the rocket without apprec- 
ieble velocity. The propulsive action, however, will remain 
the seme if the whole cartridge, that is, the propellant | 
charge end its container, is considered wholly as propellent 
but leaving the motor at a reduced effective exhaust velocity 


Al e The rocket propelled by constant thrust loses 
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only the mass equal to the propellant carried, therefore, it 
can be said to be equivalent to the “successive impulses" 
rocket if its effective exhaust velocity and its total mass of 
propellant are equal respectively to the reduced exhaust 
velocity and to the sum of the masses of all the containers of 
the "successive impulses" rocket, In other words, o is equal 
to dC and G isequito Ç’ . 

In Table 1.1 the heights for four cases have been 
calculated to illustrate the effect of the exhuast gas velocity l 
and the total number of impulses given to a rocket whose 
weight ratio, G is 0.70. It will be noticed that for 
flight in vacuo a greator hoight will be reached if a smaller 
number of impulses is employed. The lower portion of the 
Table shows the maximum height reached by an equivelent "constant 
thrust" rocket for the same four cases with the initial 
acceleration given by Eq. (lj.12c). The close agreement between l 
the maximum height reached by use of successive impulses, when 
the total numbor of impulses oxcoodo 100, and that roaohod by 
the use of constant thrust simplifies the solution of the 
problem of decressing acceleration of gravity with height, and 
enables prediction for flight with air resistance to be based 
on the results obtained for a rocket propelled by constant 
thrust (c.f. Ref. 6). These problems are considered in the 
following sections. 

Section (11) 
It is well-known that the acceleration of gravity 


decreases with tho height above tho earth's surface according 


BO 


TABLE 4.1 


. Successive impulses 


Inu, ET aa h DN] 


paso Ft. /sec. Z d y LESS Anaxo 
Ac ? per sec. Feet 
1 . 10,000 0.70 326 3 1,472,000 
2. 10,000 0.70 10 0. 092 1,686,000 
3 7,000 0.70 326 3 560,000 
4 7,000 0.70 10 0. 092 676,000 


Constant Thrust 


Has, = G EVI tapó -3)]"+ ea Leg +e 


Case Ft, /sec. Bo Emaxo 
C Č Feet 
1 10,000 0.70 32.2 1,468,000 
2 10,000 0.70 32,2 1,468,000 
3 7,000 0.70 12.9 555,000 


4 7,000 0.70 12.9 555,000 


? PRUTUTPTDÉADRDADÉDAPADRDRDRDRDRÁDÊEPITPIPEIRRRDRDIÀÁ UZ]]]]] TT EE" 


to the following relation 
Aye 
= ol. 
Ge Ga (14.33) 


At an altitude of 1000 miles the acceleration is 
only 0.6) times that at sea level. Therefore, for flights up 
to such altitudes the assumption that g is approximately 
constent is no longer valid. It was shown by Malina and Smith 
(Ref, 6) that a three-step rocket could theoretically reach such 
an altitude. Thus it is interesting to see how the decrease of 
g€ will increase the maximum height reached by the rocket, 

First the effect on powered flight in vacuo will be 
considered and then on coasting flight in vacuo, For uated 
flight the analysis is based on the assumption that the thrust 
is constent. However, the results can be applied to the case 
of propulsion by successive impulses if the total number of 
impulses, N, exceeds 100 as was justified in the previous 
section. 

The equivalent mass of gas flowing per second con- 
tinuously for the case of successive impulses is 

wn pp 
J, (4.1k) 
Assuming that the rocket starts from rest at sea 


level the equation of motion in vacuo is 
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(4.15) 


This is a non-linear differential equation which can 
not be solved by usual means. However, for all practical pur- 
poses the ratio 4 during powered flight is much smalier 
then 1, therefore, only first order terms in Ż occurring 
in the expansions need to be retained, This approximation 


linearizes the equation to the form 


ah _ 4 (4-1) He 
die do "Ta (44.16) 


g 


The solution of this equation with the initial condition that 


h=0 and Hf =o when t = O is 
Fu 
Fu E” -tu tf 
h- £] 1-05 2h Est at) -e j| EAK 
ALE e (1.17) 
em ETR ewe 
At the end of the powered flight, the time is 
pee Mg 
(4.18) 
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Therefore, the height at the end of the powered flight is 


gls) UE) 


p= Elio E Met slg) ee -Elr3) Je 
4.19) 


If the hyperbolic cosine term and the integrals are 


expanded and only first order terms iu +4 are relalned in 


consistency with the linearization of Eq. (4.15), the equation 


becomes 


y ~ - [52095 Ste Mya oT M his Megl-syesf 
+ ao (1-8) tali) 45 (18-156 ref (14.20) 


RR Le Mey Ef eli-a) togliete luct 15316] 
- apy) 


Differentiating Eq. (11.17), and substituting the 
relation of Eq. (4.18), the maximum velocity at the end of 


powered flight ia 


£lrs) 


TE val (7B dos “Jus 
| 7 Aa sah e Ms. JE: cl : IE 
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Again expending end retainiug only first order terms ln ze E 


Eq. (4.21) becomes 


la” - [Berd sm. 
É / e Laity attt 157] 


(Le22) 


-Yu AED, irofigtora aef 


It is seen that the second terms of Eq. (1.20) and 
(4.22) are the correction to be applied to Hp, and Vinay, 
to account for the variation of the acceleration of gravity. 
Since both corrections are first order approximations, they can 
be expected to epply approximately also to the case of successive 
impulses, even when the total number of impulses is less than 
100, 

The coasting height reached by the rocket due to its 
velocity at the end of powered flight can be obtained by 
equating the inorease of potential energy during coasting 


flight to the kinetic energy at the end of powered flight. Thus 


ho ti, 
hye E dh 
2 ‘may a 2 
4, (14) 


e ae 


or 
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Putting Y, , Ja (eg. ) = A which is 


coasting height obtained by assuming a constant gravitational 
acceleration of the value equal to that at the height Y. 
i.e., the height where coasting starts, then Eq. (4.23) can 


ba written 


/ 
A= (AAR oy 


MAR 


Upon expanding the second term this equation becomes, 
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This equation shows that 1f the coasting flight 
starts from sea level, and if the maximum height reached is 
about 1000 miles, the increase due to the decrease in g is 


over 25%, which is considerable, 
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Section (IT) 

When the sovnding rocket is ascending through the 
atmosphere instead of in vacuo, air resistance comes into 
play, causing the acoeleration of the rocket to be reduced, 
which decreases the maximum height reached. Since air resis- 
tance increases with the air density and with the square of 
the flight velocity, it is desirable to keep the rocket from 
ascending too rapidly through the lower layers of the atmos- 
phere where the air density is high. For this reason the op- 
timum initial acceleration will no longer be infinite as shown 
by Eq. (4.12b). For the case of constant thrust Melina and 
Smith (Ref. 6) have found thet the optimm acceleration is 
around 30ft./sece2, For a total number of impulses greater 
then 100, the difference between propulsion by successive 
impulses and by constant thrust is very small, so one may 
expect the above optimum value of initial acceleration to hold 
for both cases of propulsion. 

The actual amount of reduction in maximum hei ght 
due to air resistence can be calculated by the method of 
step-by-step integration, if fair accuracy is desired. This 
integration is carried out by using the fundamental equation 
for vertical rocket flight which, as given in the previous 


paper (Ref.6) is 
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The significance of the ratio WV. was 
0 


discussed in that paper (Ref. 6). Greater significance can, 
however, be allached Lo the various terms in the equation if it 


is transformed into the non-dimensional form 
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In Eq. (1.26) appear two types of significant 
| quantities. First, quantities, called “factors”, which ure 
constent for any given family of rockets, and second, two 
quantities celled "parameters", one of which is characteristic 
for a given family of rockets but ohanges in value along the 
flight path, and one which depends on the physical properties 
of the atmosphere. Thus there are the following factors: 
_& = ratio of initial acceleration to 


"initial acceleration factor", a motor 


o characteristic 


C = exhaust velocity in ft./sec. ~ "exhaust 
velocity factor", a motot characteristic 


A = "drageweight factor" 


Č = ratio of weight of combustibles to totel 
initial weight of the rocket ~ “loading 
factor” 


The first two factors, 1.6., the "initial secelera- 
tion factor” and the "exhuast velocity factor", determine the 
cheracteristics of the propelling unit for a given family of 
rockets while the "drag-weight factor" and the "loading factor" 
determine the physical dimensions of the rockets. The "dreg- 
weight factor” is a ratio of the drag of the rocket at sea 
level when traveling with the velocity of sound to the initial 
weight of the rocket. Since for any given family of rocket 
shapes the only terms in the factor which can be varied are the 
maximum cross-sectional area A, and the initial weight Wo, it 
is clear that if the initial weight is doubled then the cross- 
sectional area must also be doubled to keep the factor the same, 
The"loading factor needs to be discussed in some detail as it 
does not appear explicitly in Eq. (4.26). The Eg.(1,.26) is a 
differential equation of the flight path which is AO at 
every point along the flight path. The loading factor g 
comes in only when this equation is integrated and the limits 
of integration ere put in. For exemple, consider two rockets 
with identical performance factors and parameters, with the 
exception that one has a E of 0.90 and the other has a 


É of 0.50. The flight path of the two rockets will be 


«BB. 


identical up to the time that 0.50 times the initial weight of 
the rockets is used up as combustibles, At this point the 
rocket having a Z, of 0.50 will begin to dacelerata while 
the one having a č of 0.90 will continue to accelerate 
until the remaining combustibles are used up. It is thus seen 


that the value of Z controls the maximum height reached. 


The two performance parameters sre: 


q ~ physical properties of the atmosphere called 
To the "atmosphere paremeter” 

6 ~ serodynamic properties of the rocket shell 
Ge called the "form parameter" 


The "atmosphere parameter” for the earthB atmospheric - 
layer will, of course, be the seme for all rockets if standard 
conditions are assumed and its value depends only on the height 
the rocket has reached above the starting point of the flight. 
The "form parameter" is determined by the shape of the curve of. 
Cy ageinst B, This curve will be altered chiefly by the 
geometrical shape of the shell although it is also afrected by 
the change in skin friction coefficient due to the change in 
Reynold's Number. As long as the rocket belongs to a family 
that has the seme geometrical shape, which implies the same 
nose shape and the same L£/¿ ratio, that is, the ratio of 
the length of the shell to the maximum diemeter, the "form 


paremeter" can be assumed to remain constant. 
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It is thus seen that the performance curves calculated 
for a typical rocket will also hold for a whole family of rockets 
determined by the values of the "factors" and of the "parameters" 
of the typical rocket and the design of a rocket to meet certain 
prescribed requirements is greatly simplified. Furthermore, for 
a good rocket form design the variation of —> B* , the 
form parameter, at the same values of B is small. Also, the 
deviation from standard atmospheric characteristics cannot be 
very large, Then, in view of the fairly accurate but not exact 
basic avoumption of constant thrust, it is justified te use the 
seme data for these two parameters for all cases. Thus, the 
performance problem is further simplified and depends only upon 


Ga a 
the four performance factors 4 E e T AN and G , 
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CONCLUSION 


This study shows that a sounding rocket propelled 
by successive impulses can theoretically reach heights of 
much use to thoseinterested in obtaining data on the structure 
of the atmosphere and extra-terrestial phenomena if a propelling 
unit gives an exhaust velocity of 7000 ft. per second or more. 

The possibility of obtaining such exhaust velocities 
depends on two factors: first, the ability of the motor to 
transform efficiently the heat energy of the fuel into kinetic 
energy of the exhaust gases, and secondly, the amount of heat 
ener gy that een be liberated from the fuel. In an actual motor 
which burns its fuel at constant volume by igniting a powder 
charge in the combustion chamber the ratio of the chamber 
pressure to the outlet pressure drops from a maximum at the 
beginning of the expansion to zero at the end of the process. 
it is not possible to design a nozzle that will expand the 
products of combustion smoothly during the whole process. 
Therefore, the attainable efficiency must be less than that 
of a corresponding “constant pressure” motor which has a mixture 
of combustibles, e.g. gasoline and liquid oxygen, fed continu- 
ously into the combustion chamber at a constant pressure equal 
to the maximum pressure of the “constant volume” motor. 


However, very high maximum chember pressures (up to 60,000 lbs. 
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per sqein.) can be developed in a motor using constant volume 
burning, while the chamber pressure of a motor using constant 
pressure burning is limited to much lower pressures by the 
difficulty of feeding the combustibles. Therefore, the 
efficiency that can be obtained frommtors using either of 
these processes should not be very different. As to the heat 
thet can be liberated per unit mess of fuel, the present fuel, 
such as PP powder for a constant volume motor, 

is much lower them the liquid combustibles such as gasoline 
and oxygen for a constant pressure motor. 

These considerations indicate that the attainable 
exhaust velocity of a "constant volume" motor for propulsion 
by successive impulses will probably be lower then that of a 
“constant pressure” motor for supplying a continuous Lhrusbe 
This is the reason why many experimenters abandoned the 
"constant volume” motor and turned to the "constant pressure” 
motor, the so-called liquid propellant motor. Theoretically, 
this defect of the "constant volume” motor can be compensated 
if a small total number of impulses (c.f. Fig. 4.3) is used. 
However, the use of few impulses is of doubtful practical 
value because the resulting extreme accelerations will be 
hermfúl to instruments carried and will necessitate a heavier 


construction of the rocket, 
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However, even with the lower exhaust velocities of 
the "constant volume" motor it is shom by the analysis in this 
paper Lhat with the exhaust velocity of 7000 ft./300. obtainod 
experimentally by Re H. Goddard (Ref. 1) it should be possible 
to build a poder rocket capable of rising above 100,000 feet. 
Thus it seems to the author that a rocket propelled by 
successive impulses has useful possibilities end further ex- 


perimental work is justified. 
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Fig. 1.3 


heat is transferred to wall 


Pig. 1.4 Velocity and temperature distribution when wall 
temperature is 1/4 of the free stream temperature 
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Fig. 1.5 Heat balance when the wall temperature is 
1/4 of the free stream temperature 
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